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Abstract 

The NDA C and ND*A C form factors are evaluated in a full QCD sum rule 
calculation. We study the double Borel sum rule for the three point function 
of one meson one nucleon and one A c current up to order six in the operator 
product expansion. The double Borel transform is performed with respect to 
the nucleon and A c momenta, and the form factor is evaluated as a function 
of the momentum Q 2 of the heavy meson. These form factors are relevant 
to evaluate the charmonium absorption cross section by hadrons. Our results 
are compatible with constant form factors in these vertices. 
PACS numbers 14.20.Lq 12.38.Lg 14.65. Dw 



QCD predicts that at high energy density, hadronic matter will turn into a plasma of 
deconfmed quarks and gluons, the quark-gluon plasma (QGP), and there is a great deal 
of expectation in the comunity that the QGP will be observed in heavy ion collisions at 
the relativistic heavy ion collider (RHIC), which just started to operate at the Brookhaven 
National Laboratory. As a matter of fact, very recently the NA50 Collaboration reported |lj 
that the anomalous suppression of J/ip observed in Pb+Pb collisions at CERN-SPS indicated 
already the formation of QGP, since these experimental data ruled out conventional hadronic 
models for the J/ip suppression. 

However, there are calculations that reproduce the new NA50 data reasonably well up to 
the highest transverse energies based on hadronic J/ip dissociation alone. Therefore, 

while there are suggestions that the anomalous suppression may be due to the formation of 
the QGP, other more conventional mechanisms based on J /if} absorption by comovers and 
nucleons have to be still considered. 

The main ingredient in the calculations based on hadronic J /if) dissociation is the mag- 
nitude of the J/ if) absorption cross section by hadrons, which is not known experimentally. 
In refs. |2p44| this cross section was introduced as a free parameter. Therefore, in order to 
use the J /ip suppression as a signature for the formation of the QGP in heavy ion collisions 
it is important to have a better knowledge on the interactions between charmonium states 
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and co-moving hadrons. 

Various approaches have been used in evaluating the charmonium absorption cross sec- 
tion by hadrons. Some of them use meson exchange models based on hadronic effective 
lagrangians |5|-|9|]. The couplings and form factors needed in these effective lagrangians are 
not phenomenologically known, and sometimes SU(4) relations are used to estimate them. 
However, we don't have any evidence that SU(4) relations should be taken seriously into 
account. Also, in general, it is assumed a monopole form factor at the hadronic vertex J7|||, 
with another unknown parameter (the cut-off). The results obtained for the cross section 
are very sensitive to the couplings and to the form factors [0-0. As an example, in ref. [0 it 
was shown that the nJ/ip — > DD* cross section may vary by almost one order of magnitude 
depending on the value of the cut-off. Therefore, it is very important to estimate these 
couplings and form factors with a more theoretical approach. 

In this work we use the QCD sum rule (QCDSR) method (Kj based on the three-point 
function to evaluate the NDA C and ND*A C form factors and coupling constants. In a 
previous work JTT|] we have already evaluated the NDA C coupling using the QCDSR. We 
have obtained gNDA c = 6.7 ±2.1 while SU(4) relations used in ref. || give Qndk c = 14.8. Of 
course the relevance of this difference can not be underestimated since the cross section is 
proportional to the square of the coupling constant. In this work we extend the calculation 
done before [if] and study the Q 2 dependence of the form factor for an off-shell heavy meson, 
which is the situation relevant for the calculations performed in refs. ||||]. 

The three-point function associated with a ND(D*)A C vertex is constructed with the 
two baryon currents, r]\ c and t/jv, for A c and the nucleon respectively, and the meson D(D*) 
current, j^j^) which we generically call ju'- 



A M (p,p',q) = / d A xd A y(0\T{ VA Sx)jM(y)r} N m\0)e ip ' x e 



iqy 



(1) 



with the currents given by p~2"|,[13 



VAc = ^abc{ulC^ b d b )Q c 



(2) 



(3) 



(4) 



3 n = QlvV 



(5) 



where Q, u and d are the charm, up and down quark fields respectively, C is the charge 
conjugation matrix and q = p' — p. 

The general expression for the vertex function in Eq. (0) has four independent structures 
in the case of A 5 [TT] and twelve in the case of WM. In principle any of the four (twelve) 
invariant structures appearing in the Am expression can be used to calculate the form 
factor and the sum rules should yield the same result. However, each sum rule could have 
uncertainties due to the truncation in the OPE side and different contributions from the 
continuum. Therefore, depending on the Dirac structure we can obtain different results due 



2 



to the uncertainties mentioned above. The tradicional way to control these uncertainties, 
and therefore to check the reliability of the sum rule, is to evaluate the stability of the result 
function of the Borel mass. 
Recently, in ref. 
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it was pointed out that a better determination of g n NN can be 
done with the help of the 750^ structure, since this structure is independent of the effective 
models employed in the phenomenological side and it gets a smaller contribution from the 
single pole term coming from N — > N* transition. This was confirmed also in the case of the 
9nky coupling constant |Tj|. Therefore, in this work we will also employ the 750"^ structure 
to evaluate the NDA C form factor and coupling constant, and compare our results with the 
previous evaluation [IT] carried out in the ^75 structure. In the case of ND*A C vertex, we 
will study the sum rule based on the fli^p" structure since, in general, the structures with a 
large number of 7 matrices are more likely to be stable. 

The phenomenological side of the vertex function is obtained by the consideration of the 
A c and N intermediate states contribution to the matrix element in Eq.(ffl): 



Aphen), x (f + M A cl , / 



2 (_1 + Mn)_ 



+ higher resonances 



(6) 



where Aa c and \n are the couplings of the currents with the respective hadronic states and 
the meson couplings are given by 

m 2 D f D g N DA c (q 2 ) 



9s(q 



lib- 



rae q 



m 



2 ' 
D 
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9M 2 ) = m D*fo* 



9ND*kM 2 



q 
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D* 



m 
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D* 



where mp, m* D , fo and fn* are the mass and decay constant of the mesons D and D* 
respectively, and m c is the c quark mass. 

We will write a sum rule in the structure cr^j^p'^pu for gNDA c (q 2 ) and one in the structure 
$1$ f° r 9ND*A c {q 2 ) and we call F the invariant amplitude associated with these structures. 
The contribution of higher resonances and continuum in Eq. ([]) will be taken into account 
as usual in the standard form of ref. 
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In the OPE side only even dimension operators contribute to the chosen structures, 
since the dimension of Eq.([|) is four and p'p takes away two dimensions. The diagrams that 
contribute, after a double Borel transformation, up to dimension six are shown in Fig. 1. 
The gluon condensate also contributes, but it always appears with a large suppression factor 
which arises from the two-loop internal momentum integration. Therefore, its contribution is 
of little influence and will be neglected. To evaluate the perturbative contribution (Fig. la) 
we write a double dispersion relation to the invariant amplitude, F, over the virtualities p 2 
and p' 2 holding Q 2 = —q 2 fixed, and use the Cutkosky's rules [18[] to evaluate the double 
discontinuity (see ref. |T7|). After doing a double Borel transformation [17] in both variables 



P 2 = —p 2 — > M 2 and P' 2 = —p' 2 — » M' 2 , and subtracting the continuum contribution, we 

get 



F(M 2 



M' 2 " 2 



Q) 



An 2 
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-u/M' 2 -s/M 2 
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with 

p(u,s,Q 2 ) 



+ 



2m 4 Q 2 ^ 



dm 2 < m 2 



m 2 c (s - u - Q 2 ) + (s - u) 2 + Q 2 (s + u) 



A(s, w, Q 2 
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9(1 - (cos^) 2 )6 [u-Q 2 -m 2 c + 
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where 



cos 9k 



2s 



u + m 2 — m 2 — Pq(s + m 2 )/ \fs 



m*)J\(s } u } QZ) 



(11) 



with p' = (s + u + g 2 ) /(2 V«) and X(s, u, Q 2 ) = s 2 + u 2 + Q 4 - 2su + 2Q 2 s + 2Q 2 u. 

In Eq.@ F stands for the double Borel transformation of the amplitude F, and the 
subscript a refers to the diagram in Fig. la. uq and sq give the continuum thresholds for the 
baryons A c and nucleon respectively, which are, in general, taken from the mass sum rules. 

The next lowest dimension operator is the quark c mass times the quark condensate 
with dimension four (Fig. lb). Since we are neglecting the light quark masses, only terms 
proportional to m c (qq) will appear. These terms give, after the double Borel transformation 
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F(M 2 ,M' 2 ,Q 2 ) 
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The next contribution comes from the diagram involving dimension 6 operator of the 
type (qqqq) (~ (qq) 2 ) shown in Fig. lc: 



F(M 2 ,M' 2 ,Q 2 ] 



(QQ) „-ml/M' 2 



(14) 



In Eqs. (0) (pj) and QTJD the + and - sign refers to the NDA C and ND*A C vertices 
repectively. 

The Borel transformation of the phenomenological side gives 



F(M 2 ,M ,2 ,Q 2 ) 
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where the continuum contribution has already been incorporated in the OPE side, through 
the continuum thresholds sq nd uq. In Eq. fll~5|) we have defined 



Gnma c {Q' 
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(16) 



In order to obtain gND(D*)A c (Q 2 ) we identify Eq. (fl5|) with the sum of Eqs. (|j), (|12j) and 
(0). We obtain: 



e M */ M2 e M V M ' 2 m c (Q 2 + ml) 
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du [° dse~ u ' M ' 2 e- s l M2 (p(u,s,Q 2 ) 

(18) 



For Aa c and Ajv we use the expressions obtained from the respective mass sum rules for 
uc 



A c [13] and for the nucleon 12,19 
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\X N \ 2 = e M ^^E 2 + l(qq) 2 ^ , (20) 

where E 2 = 1 — e _S0//M Ai(l + so/Mf^ + Sq/(2M^-)) accounts for the continuum contribution. 
For consistency we have also neglected the contribution of the gluon condensate in the mass 
sum rules, since it is of little influence. 

In Eqs. (HD and (||) M 2 M and M' M represent the Borel masses in the two-point function 
of the nucleon and A c respectively. Comparing Eqs. (|T7|), QXD|) and (^) we can see that the 
exponentials multiplying Eq. ([[7]) disappear if we choose 

2M 2 M = M 2 and 2M' M 2 = M' 2 . (21) 

Indeed, this way of relating the Borel parameters in the two- and three-point functions 
is a crucial ingredient for the incorporation of heavy quark symmetries, and leads to a 
considerable reduction of the sensitivity to input parameters, such as continuum thresholds 
so and uq, and to radiative corrections [|l~4"l , p0l . 



The parameter values used in all calculations are m c = 1.5 GeV, = 1.87 GeV, 
m D * = 2.01 GeV,MN = 938 MeV, M Ac = 2.285 GeV, f D = 170 MeV f D * = 240 MeV 0, 
(qq) = — (0.23) 3 GeV 3 . We parametrize the continuum thresholds as 

s = (M N + A s ) 2 , (22) 
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and 



u = (M Ac + A u ) 2 . (23) 

The values of uq and so are extracted from the two-point function sum rules for M\ c and 
in Eqs. ( |19D and (|20"D and the respective sum rules in the 1 structure given in refs. []T2| , |T3 , 1 9 



We found a good stability for M^ c and M N , being able to reproduce the experimental values 
for the masses in the Borel mass M 2 ^ ~ 1 GeV 2 and M' M 2 ~ 6GeV 2 , with A s = 0.7 GeV 
and A u = 0.6 GeV, which are the values that we are going to use in the calculations. 

To allow for different values of M 2 and M' 2 we take them proportional to the respective 
baryon masses. In this way we study the sum rule as a function of M 2 at a fixed ratio 

Ml = Mk (24 ) 
M' 2 M 2 c ' [ ' 

In Fig. 2 we show the behavior of the perturbative, quark condensate and four quark 
condensate contributions to the form factor gNDhXQ 2 ) at Q 2 = GeV 2 (since in this case 
the cut in the t channel starts at t ~ m 2 and thus the Euclidian region stretches up to that 
threshold) as a function of the Borel mass M 2 . We observe that the different contributions 
add to give a very stable result as a function of the Borel mass. Since the Borel masses in the 
two- and three-point functions are related by Eq. (^l|) and since ~ 1 GeV 2 , to study the 
Q 2 dependence of the form factor we fix M 2 = 2.5 GeV 2 where the perturbative contribution 
is the dominant one. The behaviour of the curve for other Q 2 values is similar, however, 
for Q 2 > 1.5 GeV 2 the perturbative contribution is no longer the dominant one. Therefore, 
in Fig. 3 we show the Q 2 dependence of the form factor in the region < Q 2 < 1.5 GeV 2 
(dots) where we believe we can trust the QCDSR results. We see that the form factor is 
practically constant in this region, showing a tiny decrease. We can fit the QCDSR results 
with a monopole form, as can be seen by the solid line in Fig. 3, and we get gNDhXQ 2 ) = 
843/(102+Q 2 ) which corresponds to a cut-off of order of 10 GeV, much bigger than the values 
used in ref. ||. If we vary the value of the Borel mass used to extract the Q 2 dependence 
of the form factor we can even get a result that shows a tiny increase in the considered 
region. Therefore, in view of the uncertainties involved in the approach, we can say that 
the QCDSR give a constant value for the form factor. Considering 20% variation in the 
continuum thresholds and fo varying in the interval fo = 170 ± 10 MeV [[H]] we get: 



dNDKciQ 2 ) = 7.9 ± 0.9 , (25) 



in agreement with our previous estimate of the coupling constant Jy 

The same analysis can be done for gND*A c (Q 2 ) and in Fig. 4 we show the perturbative, 
quark condensate and four quark condensate contributions to the form factor gND*A c (Q 2 ) 
at Q 2 = GeV 2 as a function of the Borel mass M 2 . Since the OPE sides of both form 
factors differ only by a sign the aspect of Fig. 2 is very similar to Fig. 4. In Fig. 5 we show 
the Q 2 dependence of the form factor gND*A c (Q 2 ), extracted at M 2 = 2.5 GeV, in the region 
< Q 2 < 1.5 GeV 2 (dots). Again we observe that the form factor is practically constant 
in this region, showing a tiny decrease. Fitting the QCDSR results with a monopole form 
(solid line in Fig. 5) we get gND*A c (Q 2 ) = —204/(26.2 + Q 2 ) which corresponds to a cut-off 
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of order of 5 GeV, still much bigger than the values used in ref . . Despite the fact that 
for this particular choice of M 2 , s and Uq we get a form factor with a smaller cut-off than 
obtained for gNDA c (Q 2 ), varying M 2 and the continuun thresholds we still can get results 
that show a tiny increase in the considered region. Therefore, also in this case, due to the 
uncertainties involved in the approach, we conclude that the QCDSR give a constant value 
for the form factor. Considering 20% variation in the continuum thresholds and fn* varying 
in the interval f D . = 240 ± 20 MeV |2J we get: 

9nd*k c {Q 2 ) = -7.5 ±1.1 , (26) 

It is important to mention that constant form factors in hadronic vertices have already 
been found in the context of QCDSR. This is the case of the goD*n form factor. This form 
factor has been studied in the QCDSR approach for an off-shell D* pl]-p^1, and for an 
off-shell pion [23]. As a matter of fact, in refs. |2"T|-|2"3"|| the authors analyze the semileptonic 
/+(£) form factor defined by 

(ir(p v )\uj^c\D(p D )) = (p ff +p D )J+(t) + (pD-p*)pf-(t) , (27) 

where t = q 2 = (po — Pn) 2 - Since the vector current = u^^c has the same quantum 
numbers as the vector meson D*, the same sum rules studied in refs. [PTH 23H can be used 
to study the hadronic form factor gDD*w{t), for an off-shell D* meson. It is only the phe- 
nomenological side of the sum rule that has to be modified to allow for a coupling between 
the vector current and the vector meson: 

(tt^I^ID^)) = (0\V,\n(- P7T )D(p D )) 

= (0K\D*( q ))-^-^(D*( q )n( P7T )\D(p D )) . (28) 

The gDD*n(t) form factor is defined by the strong amplitude 

(D*(q,e)n(p T )\D(p D )) = gnn^pn + p T ) a e a . (29) 

Therefore, using the the vacuum to vector meson transition amplitude defined in terms of 
the vector meson decay constant fr,*: 

(D*{q,e)\V^)=m D *f D ,el- (30) 

we can rewrite Eq. fl2jf) as 



{<Pir)K\D{p D )) = —— 2 (p n + p D )^ - . (31) 

i m D * \ ra D * j 

Comparing Eqs. (|27| ) and (|3~ID we imediately see that goD*^) and /+(£) are related by 

f + (t) = m ^9DD^(t) (32) 
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The authors of refs. jn],[22| claim that the QCDSR results for f+(t) can be well fitted 
by a monopole form with a pole mass m po i = mp*. According with Eq. (|32"D this implies a 
constant gDD*n(t). More recently the authors of [23|, tried to fit their QCDSR results with 
a double pole parametrization of the type 

= (l-t/ml.)(l- at/ml.) ' (33) 

They have obtained a = 0.01 and, therefore, they have concluded that a is consistent 

with zero, which means the complete dominance of the vector meson pole to / + with a 
consequent prediction of a constant hadronic form factor in the DD*n vertex. 

It is interesting to notice that, due to the uncertainties in the value of a given above, 
one can not even rule out a DD*ir form factor that slightly growns in the Euclidian region, 
as mentioned above in the case of ND(D*)A C . 

A very different result for the form factor in the DD*n vertex was obtained in the 
case of an off-shell pion EM. This form factor is important in the processes described in 



ref. ||25|| ) which could contribute to explain the enhacement in the production of dileptons of 
intermediate masses, also observed in the NA50 experiment. In the case that the pion is off- 
shell, the form factor shows a very pronounced Q 2 dependence, with Q 2 being the squared 
momentum of the off-shell pion in the Euclidian region. This result might be suggesting that 
the size of a hadronic vertex depends on which particle is off-shell. If the off-shell particle 
is light, then the vertex is not point like. However, if the off-shell particle is heavy, then the 
vertex is point like with a consequent constant form factor. Our result for the ND(D*)A C 
form factors corroborate this hypothesis. 

In conclusion, in this work we have calculated the form factors for the hadronic vertices 
NDA C and ND*A C using QCD sum rules. These form factors are important to evaluate 
the charmonium dissociation cross section by hadrons, in the framework of meson exchange 
models based on effective lagrangians. In the construction of the sum rules we have per- 
formed a double Borel transformation with respect to the nucleon and A c momenta, and we 
have evaluated the form factors as a function of the heavy meson momentum Q 2 . We have 
studied the sum rules in the structures cr^j^p'^Pu for gND\ c (Q 2 ) and ft^rfj) for gND*K c {Q 2 )- 
In the studied Q 2 region, our results are compatible with constant form factors in these 
vertices. Considering 20% of variation in the continuum thresholds, and around 10% of 
variation in the meson decay constants we got: 

2vda c (Q 2 ) = 7.9±0.9 

g N D*A c (Q 2 ) = -7.5 ±1.1 . (34) 

It is important to stress that constant form factors where also found, in the framework 
of QCDSR, in the vertex D*Dn with an off shell D*. Off course for very large values of Q 2 
the form factors should go to zero, but, according to our results, this means that the values 
of the cut-offs in these form factors should be much larger than 2 GeV, which is the typical 
value used in refs. ]7|-[| . 
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FIG. 1. Diagrams that contribute to the OPE side of gND(D*)A c (Q 2 )- The cross stands for a 
charm mass insertion. 
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FIG. 2. Borel dependence of the perturbative (dashed line), quark condensate (dotted line) 
and four quark condensate (dot-dashed line) contributions to the NDA C form factor (solid line) at 
Q 2 = 0. 
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FIG. 3. Momentum dependence of the NDA C form factor (dots). The solid line give the 
parametrization of the QCDSR results with a monopole form. 
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FIG. 4. Borel dependence of the perturbative (dashed line), quark condensate (dotted line) 
and four quark condensate (dot-dashed line) contributions to the ND*A C form factor (solid line) 
at Q 2 = 0. 
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FIG. 5. Momentum dependence of the ND*A C form factor (dots). The solid line give the 
parametrization of the QCDSR results with a monopole form. 
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